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We study the scaling of two-dimensional crack roughness using large scale beam lattice systems. 
Our results indicate that the crack roughness obtained using beam lattice systems does not exhibit 
anomalous scaling in sharp contrast to the simulation results obtained using scalar fuse lattices. 
The local and global roughness exponents {(loc and (, respectively) are equal to each other, and 
the two-dimensional crack roughness exponent is estimated to be (loc = C = 0-64 ± 0.02. Removal 
of overhangs (jumps) in the crack profiles eliminates even the minute differences between the local 
and global roughness exponents. Furthermore, removing these jumps in the crack profile completely 
eliminates the multiscaling observed in other studies. We find that the probability density distribu- 
tion p{Ah{£)) of the height differences Ah{£) = [h{x + £) — h{x)] of the crack profile obtained after 
removing the jumps in the profiles follows a Gaussian distribution even for small window sizes {£). 



I. INTRODUCTION 



Understanding the scaling properties of fracture sur- 
faces still represents an unsolved problem despite two 
decades of intense research activities [l|, 0] ■ Experiments 
on several materials under different loading conditions 
have shown that the fracture surface is self-afhne [3] and 
can be characterized by a roughness exponent C- Ex- 
eriments on several materials including metals [3] , glass 
rocks (6| and ceramics [3], have shown a universal 
out of plane roughness exponent of ^ ~ 0.8 for three- 
dimensional fracture surfaces irrespective of the material 
studied, as revieved in Ref. Q- Recent experimental 
evidence shows, however, that the picture is more com- 
plicated: such scaling is valid at small and intermediate 
scales in the so-called fracture process zone (FPZ), while 
at large scales one observes a new regime with C ~ 0.4 
attributed to crack line depinning [3, \1Q, [ll| . Recent nu- 
merical investigations of the random fuse model indicate 
that the local roughness exponent in two dimensions does 
not depend on the size of the FPZ, but only on the fact 
that a FPZ is present or not [12| . 

In addition, it is by no means a priori clear that "sim- 
ple" self-affinity is sufficient to describe the experiments. 
It has been argued that fracture surfaces may exhibit 
anomalous scaling [l3|: the global exponent describing 
the scaling of the crack width with the sample size is 
larger than the local exponent measured on a single sam- 
ple [ij, \lM- ^^ t^i^ sense, it is necessary to introduce 
two roughness exponents a global exponent (C) and a lo- 
cal exponent (Ooc) to define the roughness of fracture 
surfaces. Anomalous scaling is noted in the numerical 



simulations as well [la . [r7| : however, its origin is not 
clear yet although in experiments it is conjectured to be 
an artifact of initial transient regime as the fracture front 
moves away from the initial notch [Il|. Recent studies 
[Wi . [ig | also suggest that the origin of anomalous scaling 
in numerical simulations in two dimensions may also be 
due to the existence of overhangs (jumps) in the crack 
profile, originating from crack branching. As a further 
complication, there is an on-going debate whether frac- 
ture surfaces exhibit multi-affine scaling [1^, [l9|, [20|, |2l| , 
implying that one would observe a whole family of rough- 
ness exponents ^q depending on which statistical moment 
q of the correlation function is measured. 

The theoretical understanding of the origin and uni- 
versality of crack surface roughness is often investigated 
by discrete lattice (fuse, central- force, and beam) models 
y. In these models the elastic medium is described by 
a network of discrete elements such as fuses, springs and 
beams with random failure thresholds. In the simplest 
approximation of a scalar displacement, one recovers the 
random fuse model (RFM) where a lattice of fuses with 
random threshold are subject to an increasing external 
voltage [22|. Using two-dimensional RFM, the estimated 
crack surface roughness exponents are: ( = 0.7 ± 0.07 
[H, (i^^ = 2/3 in, and C = 0.74 ± 0.02 [25|. Recently, 
using large system sizes (up to L ~ 1024) with exten- 
sive sample averaging, we found that the crack roughness 
exhibits anomalous scaling [2(a] . In particular, the lo- 
cal and global roughness exponents estimated using two 
different lattice topologies are: Qoc — 0.72 ± 0.02 and 
C — 0.84 ±0.03. The reasons behind the origin of anoma- 
lous scaling in numerical simulations are not yet clear, 
although the existence of overhangs in the crack pro- 



file is expected to have contributed to anomalous scaling 
[la . [l9| . In comparison, the roughness exponents ob- 
tained from quasi two-dimensional experiments, mainly 
obtained for paper samples, indicate a roug hness expo- 
nent in the range C =i 0.6 - 0.7 [H, M, Hi HI, M, iS, 
but occasionally significantly higher values have also been 
reported |3l|. It is not known at this time whether this 
variation in C values is a reflection of practical difficulties 
in experimentally measuring C - in paper it is difhcult to 
have a scaling range spanning over several decades since 
the structure becomes three-dimensional at small scales 
(0.1 mm) and at the millimeter range the fiber length 
interferes - or that the roughness exponent is not really 
universal but depends on material parameters and the 
anisotropy of the medium. 

Despite this reasonable agreement between the numer- 
ical results obtained using two-dimensional RFM and the 
above quasi-two-dimensional experimental results, a lin- 
gering question is whether scalar representation of the 
elastic medium using random fuse models is an adequate 
representation of fracture. Moreover, it is an intriguing 
question whether the same roughness exponents as those 
obtained using RFM will be obtained using more com- 
plex random threshold central-force (spring) and beam 
models. Using central-force models, a roughness expo- 
nent of 0.65 ± 0.07 is obtained in Ref. [33 • A more 
recent study with a range of disorder strengths in ran- 
dom thresholds estimated the local and global roughness 
exponents to be in the range 0.5 — 0.67 and 0.6 — 0.85, 
respectively [3 • Using two-dimensional beam simula- 
tions, Ref. [33] estimated a typical value for the global 
roughness exponent to be 0.86; however, the roughness 
exponent is argued to be disorder dependent. Similar 
disorder dependent roughness exponents have also been 
reported recently for the RFM [34]. These results were 
obtained in the square lattice and were attributed to a 
lattice effect at low disorder. They indeed disappear for 
triangular lattice, where the r oug hness exponent is the 
same independent on disorder [l2l |. 

The questions we address in this article are the fol- 
lowing: (i) whether anomalous scaling is present in two- 
dimensional fracture simulations using beam lattice sys- 
tems, and (ii) whether roughness exponents so measured 
using beam lattice systems are in agreement with those 
obtained using simplified scalar RFM models. Along the 
way, we also address why scalar RFM models have been 
successful in representing fracture in a disordered elas- 
tic medium. Recent studies [l^, [20, [20| have shown that 
multi-affinity of fracture surfaces is an artifact of over- 
hangs (or jumps) in the crack profile. Here, we further 
investigate the influence of these overhangs in the crack 
profiles on crack roughness exponents. This article has 
three further Sections: in the next one, we describe the 
beam model used. In Section III, we present the numeri- 
cal results obtained with it. Finally, Section IV concludes 
the paper. 



II. BEAM MODEL 

The random thresholds beam model (RBM) we con- 
sider in this study is a two-dimensional triangle lattice 
system of size L x L. Unlike the scalar RFM model, 
the vectorial RBM has three degrees of freedom (x- 
translation u, y-translation v, and a rotation 9 about z 
axis) at each of the lattice nodes (sites), and each of the 
bonds (beams) in the lattice connects two nearest neigh- 
bour nodes. We assume that the beams are connected 
rigidly at each of the nodes such that the angle between 
any two beams connected at a node remains unaltered 
during the deformation process. These nodal displace- 
ments and rotations introduce conjugate forces and bend- 
ing moments in the beam members. Using Timoshenko 
beam theory [33|, which includes shear deformations of 
the beam cross-section in addition to the usual axial de- 
formation of cross-sections, the Zoca/ stiffness matrix for a 
beam element that relates the local nodal displacements 
and rotations to local nodal forces and bending moments 
in the beam's local coordinate system is given by 
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where E is the Young's modulus, G is the shear modulus, 
A is the beam cross-sectional area, / is the moment of in- 
ertia of beam cross-section, £& is the length of the beam, 
and a = l?fl> is the shear correction factor, which de- 

notes the ratio of bending stiffness to the shear stiffness. 
If shear deformation of beam cross-section is negligible, 
then a ~ and the Timoshenko's beam theory reduces 
to Euler-BernouUi beam theory. Equation [1] presents a 
relation between local nodal displacements and rotations 
di = {uii,vii,dii,uij,vij,9ij)'^ and local forces and mo- 
ments Fi = {Fii,Vii,Mu,Fij,Vij,Mi-i)'^. In this setting, 
the subscript I refers to local quantities, the superscript 
T represents transpose of a vector or a matrix, i and j 
refer to i-th and j-th nodes of the beam, and F, V, and 
M refer to axial force, shear force, and bending moments 
respectively. 

Computing the equilibrium of the lattice system is 
achieved by first transforming these local quantities (d^ 
and Fi) into global quantities d = {ui,Vi,9i,Uj,Vj,9j)'^ 
and F — {Fi,Vi,Mi,Fj,Vj,Mj)'^ through a coordinate 
transformation T such that d^ = Td, F^ = TF, and 
K = T"^K;oca;T, and then satisfying equilibrium equa- 
tions at each node such that 
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where S<ij> implies that the summation is carried over 
all the intact bonds < ij > joining at node i. In the 



above discussion, the transformation matrix T is given 
by 
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and c = cos(/3), s = sin(/3) refer to the direction cosines 
of the beam with (3 representing the angle between the 
beam axis and the x-direction. 

In this RBM, we start with a fully intact lattice sys- 
tem with beams having unit length, unit square cross- 
section and Young's modulus E — I. This results in 
a unit axial stiffness {EA/ib = 1) and bending stiffness 
{12EI/£l = 1) for each of the beams in the lattice system. 
Since the beam can deform in two independent deforma- 
tion modes (axial and bending), we assume randomly 
distributed bond breaking axial and bending thresholds, 
ta and ift, based on thresholds probability distributions, 
Pa{ta) and Pb{tb) respectively. The failure criterion for a 
beam is defined through an axial force and bending mo- 
ment interaction equation (similar to von-Mises criterion 
in metal plasticity) given by 



max(|Mi|,|Mj| 



tb 



1 



(5) 



The beam breaks irreversibly, whenever the failure crite- 
rion r > 1. Periodic boundary conditions are imposed in 
the horizontal direction and a constant unit displacement 
difference is applied between the top and the bottom of 
lattice system. 

Numerically, a unit displacement, A = 1, is applied at 
the top of the lattice system and the equilibrium equa- 
tions (Eq. [21) are solved to determine the force in each of 
the springs. Solution of Eq. [5] results in global displace- 
ments and rotations d, using which the local displace- 
ments df = Td and the local forces F^ = K.iocai<ie are 
computed for each of the intact beams. Subsequently, 
for each bond k with nodes i and j, the quantities 
max(|Mi|,|Mj 



flfe 
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and bl- 



are evaluated, and 



the bond fee having the smallest value, 



rk 



-bk + y^b'i + 4afc 
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(6) 



is irreversibly removed (When ak — 0, then rk = —)■ 
The forces are redistributed instantaneously after a bond 
is broken implying that the stress relaxation in the lattice 
system is much faster than the breaking of a bond. Each 
time a bond is broken, it is necessary to re-equilibrate 
the lattice system in order to determine the subsequent 
breaking of a bond. The process of breaking of a bond, 
one at a time, is repeated until the lattice system falls 



apart. For the RBM, we consider a uniform probability 
distribution in [0, 1] for both axial and bending thresholds 
disorders. 

Numerical simulations of fracture using lattice net- 
works have often been limited to smaller system sizes 
due to the high computational cost associated with solv- 
ing a new large set of linear equations every time a 
new lattice bond is broken. In this work, we use the 
multiple-rank sparse Cholesky factorization downdating 
algorithm developed by the authors for simulating frac- 
ture using discrete beam lattice systems [3a, |33!- For 
two-dimensional systems, this low-rank Cholesky factor 
downdating algorithm is significantly faster than com- 
peting preconditioned conjugate-gradient based iterative 
solvers. Using this numerical algorithm, we were able 
to investigate fracture in larger lattice systems (e.g., 
L = 320 in 2D) than those investigated in previous stud- 
ies. The lattice system sizes considered in this work are 
L = {32, 64, 128, 256, 320} with large numbers of sample 
configurations, Nconfzg = {2000,4000,640,400,200} re- 
spectively, in order to reduce the statistical error in the 
numerical results. 



III. CRACK ROUGHNESS 

Once the sample has failed, we identify the final crack, 
which typically displays dangling ends (see Fig. [T]). We 
remove them and obtain a single valued crack line hx, 
where the values oi x € [0,L]. For self-affine cracks, the 
local width, w{l) = {J2^{h, - (1/0 Ejf hxW^^, where 
the sums are restricted to regions of length I and the 
average is over different realizations, scales as w{l) ~ K 
for I -C L and saturates to a value W = 'w{L) ^ L'' 
corresponding to the global width. The power spectrum 
S{k) = {hkh_k)/L, where hk = X^a; ^a; expi(27rxfc/i), 
decays as Sjk) ^ fc~(^'»+^^. When anomalous scaling 
is present llSl . \14 . Il5| , the exponent describing the sys- 
tem size dependence of the surface differs from the local 
exponent measured for a fixed system size L. In par- 
ticular, the local width scales as w{t) ^ (Cioc j^Q-Qioa ^ gQ 
that the global roughness W scales as L^ with Q > Cioc- 
Consequently, the power spectrum scales as iS'(fc) ^ 

^-(2C,„c + l)^2(C-C(oc), 

Recently, Bouchbinder et al. [2l| have suggested that 
the crack line h{x) is not self-affine; instead, it exhibits 
a much complicated multi-affine (or multiscaling) struc- 
ture. This implies a non-constant scaling exponent Cg for 
the q-th order correlation fimction Cq{i) = {\h{x + i) — 

/i(a;)|9)i/'? ~ e'^" iglil. This would imply the breakdown of 

— ' I — ~\ I — 1 1 — ^ 

self-affinity. Recent studies [181, |20|, |29|| have shown that 

multi-affinity of fracture surfaces seems to disappear on 
large enough scales £, however. 

As shown in Fig. ^ removal of jumps from an initially 
periodic crack profile h(x) makes the resulting crack pro- 
file hisip{x) nonperiodic, where the subscript NP refers 
to nonperiodicity of the profiles. A direct evaluation of 
the roughness exponent using these nonperiodic profiles 
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FIG. 1: (Color online) A typical crack in a system of size LxL 
with L — 256. This crack, identified as (a) in the figure, has 
dangling ends, which are removed to obtain a single valued 
crack profile h(x), identified as (b) in the figure. However, 
this final crack h{x) possesses finite jumps that arise due to 
solid-on-solid projection of fracture surfaces. 



can be made. However, such an evaluation of roughness 
exhibits finite size effects for window sizes £ > L/2 due 
to nonperiodicity. Alternatively, the roughness of these 
resulting profiles can be evaluated by first subtracting 



a linear profile hunix) = 



INP 



(0) 



{hMp{L)-hNp(0)) . 



from the nonpcriodic profile hNp{x), and then evaluat- 
ing the roughness of the resulting periodic profile hp{x). 



show the crack widths calculated based on original crack 
profiles with jumps and those obtained from crack pro- 
files without the jumps. The jumps in the profiles appear 
to result in slightly different local and global roughness 
exponents {C,ioc and Ci respectively), although the expo- 
nents are within error bars. However, removing these 
jumps in the crack profiles leads to a single roughness 
exponent of Qioc = C = 0-65 suggesting that anomalous 
scaling is an artifact of jumps in the crack profiles, at 
least for fracture simulations based on beam lattice sys- 
tems. We have also investigated the power spectra S{k) 
of the crack profiles with and without the jumps in the 
crack profiles (see Fig. [3j)). Collapse of the power spec- 
tra for different system sizes can be observed for both 
the sets of crack profiles and the roughness exponents 
(Ooc = C = 0.65) obtained using the power law fits to the 
data are consistent with those presented in Fig. [3^. 

The self-affine property of the crack profiles also im- 
plies that the probability density distribution p{/S.h{i)) 
of the height differences /S.h{t) = [h{x ~^ t) — h{x)] of the 
crack profile follows the relation 



v{^h{i)) r. {Ah\e))-'/' f 



Ah{£) 



(A/i2(^))i/: 



(7) 



Noting that periodicity in crack profiles is analogous 
to return-to-origin excursions arising in stochastic pro- 
cesses, we propose the following ansatz for the local width 
(A/i2(£))i/2 in height differences Ah{£) 



{Ah\£)y/' =. (A/,2(^/2))i/2^(_i_ 
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(8) 
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FIG. 2: (Color online) Figure shows a typical single valued 
crack profile h{x) with jumps based on solid-on-solid projec- 
tion scheme (identified as (a)). Removing the jumps in the 
crack profile h{x) makes it a non-periodic profile (identified 
as (b)). Subtracting a linear profile from this non-periodic 
profile results in a periodic profile (identified as (c)). 



with (A/i2(L/2))i/2 = lC, For periodic profiles, the func- 
tion (j) I "ryn ) is symmetric about £ — L/2 and is con- 
strained such that (J) ( jj^ 1 = at £ = and £ ^ L, and 

(j) ( jj^ j = 1 dX £ — L/2. Based on these conditions, a 
scaling ansatz of the form 



(A/l2(^))l/2 



(A/t2(i/2))i/2 



1/Ci, 



{t-L/2f 
(L/2)2 



(9) 



similar to stochastic excursions or bridges can be pro- 
posed for (Aft,2(^))i/2^ which implies a functional form 



( £ \ \ ( {£-L/2) \ 



(10) 



for (f) I 142 ) that is satisfied to a good approximation 
by our numerical results. This scaling ansatz results in 
anomalous scaling when (lioc ^ C- Upon further simplifi- 
cation, Eq. [10] results in 



Figure [3^ presents the scaling of crack width w{£) with 
window size £. The inset and the main figure respectively 
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FIG. 3: (Color online) (a) Scaling of crack width w{(.) with 
window size £. The inset shows crack widths calculated based 
on original crack profiles with jumps. The local and global 
roughness exponents {(^loc and (, respectively) appear to be 
within error bars although the curves do not collapse. In 
the main figure, crack widths are calculated based on crack 
profiles obtained after removing the jumps in the profiles. Re- 
moving the jumps in the profiles appears to eliminate the dis- 
crepancy between local and global roughness exponents, (b) 
Scaling of power spectrum S{k). Inset shows the scaling of 
S{k) for original crack profiles whereas the main figure shows 
the scaling of S{k) for crack profiles without jumps. Collapse 
of the power spectrum can be observed, and the exponents 
are consistent with those shown in (a). 



which along with (A/i2(i/2))i/2 = ^C and Eq. [8] shows 
how anomalous scaling arises and how local and global 
roughness exponents Qoc and C can be computed based 
on numerical results. 

Figure m presents the scaling of (A/i^(^))^/^ based on 
the above ansatz (Eqs. [5] and [TU|) . The collapse of 
the (A/i|,(^))i/V(A/i|,(i/2))i/2 data for different sys- 
tem sizes L and window sizes £ onto a scaling form given 
by Eq. [TO] with (loc = 0.64 can be clearly seen in Fig. 
mja). In addition, the collapse of the data presented in 
Fig. HJc) for {Ah%{e)y/y(Ah%{L/2)y/'' provides quite 



concrete evidence that multi-affine scaling of fracture sur- 
faces, similar to what is observed in Ref. [21|, is an arti- 
fact of jumps in the crack profile that are formed due to 
the solid-on-solid approximation used in extracting the 
crack profiles. The results in Fig. [IJc) clearly demon- 
strate that the removal of these jumps in the crack pro- 
files completely eliminates this apparent multiscaling of 
fracture surfaces. This is also evident through the scaling 
of {Ah%{L/2)y/'} presented in Fig. gtb). The slopes of 
the data for moments g = 1 to 6 of Ahp{L/2) are iden- 
tical. It should also be noted that (A/i|,(L/2))i/'? - L? 
with C — 0.64 and is identical to the local roughness ex- 
ponent as obtained from Figs. Ufa) and (c). This further 
indicates that anomalous scaling of crack profiles is not 
present in fracture simulations obtained using the beam 
lattice systems. As already noted in Figs. [5^-b, the dif- 
ference between the local and global roughness exponents 
obtained using the original crack profiles is so small that 
it already negates the existence of anomalous scaling of 
crack roughness using beam lattice simulations. Remov- 
ing the jumps caused by overhangs in the crack profile 
further reduces even this minute difference in local and 
global roughness exponents thereby eliminating anoma- 
lous scaling of crack profiles. 

In the following, we investigate the probability den- 
sity »(A/i(^)) of height differences Ah{e). In Refs. 
[i3,[20,[23|, the p{Ah{£)) distribution is shown to follow a 
Gaussian distribution above a cutoff length scale and the 
deviations away from Gaussian distribution in the tails 
of the distribution have been attributed to finite jumps 
in the crack profiles. A self-affine scaling of p{Ah{£)) as 
given by Eq. [7] implies that the cumulative distribution 
P{Ah{£)) scales as P{Ah{£)) - P{Ah{£:)/ {Ah^{£))^''^). 
Figuretnija) presents the raw data of cumulative probabil- 
ity distributions P{Ah{£)) of the height differences Ah{£) 
on a normal or Gaussian paper for bin sizes £ <ti L. As 
observed in Refs. [l3, [20)[23, Fig. [5fa) shows large de- 
viations away from Gaussian distribution for these small 
bin sizes. However, for moderate £, the distribution is 
Gaussian with deviations in the tails of the distribution 
beyond the 3cr = 3(A/i^(^))^/^ limit (data not shown 
in Figure). Removing the jumps in the crack profiles 
however collapses the P{Ahp{£)) distributions onto a 
straight line (see Fig. Eljb)) indicating the adequacy of 
Gaussian distrbution even for small £. Indeed, Fig. [^b) 
shows the collapse of the P{Ahp{£)) data for systems of 
sizes L — 128 and L = 256 with a variety of bin sizes 
2 < £ < L/2. Removing the jumps in the profiles not 
only turns the P{Ahp{£)) distributions Gaussian even 
for small window sizes £ but also extends the validity of 
P{Ahp(£)) Gaussian distribution for moderate bin sizes 
to a 4cr = 4(A/i|,(^))i/2 (99.993% confidence) limit. 



IV. DISCUSSION 

In summary, the analysis and results presented in this 
paper indicate that crack profiles obtained in fracture 
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FIG. 5: (Color online) Normal paper plots of cumulative 
probability distributions P{Ah{£)) of the height differences 
Ah{£) — [h{x + £) — h{x)] of the crack profile h{x) for var- 
ious bin sizes £. $~^ denotes inverse Gaussian. The col- 
lapse of the profiles onto a straight line with unit slope would 
indicate that a Gaussian distribution is adequate to repre- 
sent P{Ah{£)). (a) P{Ah{e)) distributions for L = 256 
and £ <^ L. Large deviation from Gaussian profiles is ob- 
served for these window sizes, (b) Removing the jumps in 
the profiles however collapses the P{Ahp{i)) distributions 
onto a straight line indicating the adequacy of Gaussian dis- 
trbution even for small window sizes £. The data is for 
L = 128 with £ = (2,4,8,16,32,48,64), and L = 256 with 
£= (2,4,8,16,32,64,96,128). 
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FIG. 4; (Color online) Removing the jumps from the crack 
profiles eliminates anomalous scaling. Furthermore, multi- 
scaling of crack profiles is also due to the jumps in the crack 
profiles, (a) Scaling of {Ah%{£))^''^ with window size £ shows 
collapse of data onto Eq. [10] with (^loc = 0.64. (b) Scaling 
of (A^p(Z//2))^" with system size L with a scaling exponent 
of C = 0.64. (c) Scaling of (A/i|,(^))^''« with window size £ 
shows that the apparent multiscaling of crack profiles is due 
to jumps in the profiles. 



simulations using the beam lattice systems do not ex- 
hibit anomalous scahng of roughness in contrast with 
those obtained using the fuse lattice systems. The local 
roughness exponent Qoc = 0.65 obtained using the beam 
lattice simulations is in good agreement with the ear- 
lier results obtained using the central-force spring mod- 
els {Cioc = 0.65) [34I and is rather close to that of the 
fuse models C/oc = 0.72 [2^1 ■ Notice, however, that re- 
cently it was shown that for the random fuse model on 
the square lattice at low disorder the roughness exponent 



is larger [34|. This result was attributed to bias due to 
the lattice topology. In fact, for the trangular lattice this 
effect is not seen and the roughness exponent is universal 
|12 | . Lattice could also explain the difference between the 
present result for the beam model (Qoc — 0.65) and pre- 
vious disorder dependent results obtained for the square 
lattice [33|- 

The agreement in the local roughness exponent be- 
tween the beam, spring, and possibly fuse models is in- 
teresting to note because of these model's dissimilarity 
in representing deformation of an elastic medium. This 
seems to imply that the anisotropy in the stress redis- 
tribution in tensorial models is irrelevant for the rough- 
ness, at least in two dimensions. We even obtain the 
same local roughness exponent by considering a simpli- 
fied RBM in which failure events form a connected crack 
thereby excluding damage nucleation in the bulk. This 
simplified beam model is similar to the simplified ran- 
dom thresholds fuse model (RFM) considered in the Rcf. 
|la |. In this model, after breaking the weakest beam, 
successive failure events are only allowed on beams that 
are connected to the crack. Otherwise, the rules of this 
simplified model strictly follow those of the usual RBM. 
Consequently, this model tracks only the connected crack 
along with its dangling ends in a disordered medium, and 
hence forms the most simplified model to study the effect 
of disorder on crack roughness. 

As noted in Ref. [18| for the fuse models, this simplified 
beam model exhibits the same local roughness and power 
spectra characteristics as that of conventional RBM (see 
Fig. [S^) . This implies that one can expect to obtain same 
local roughness exponent as long as there exists a frac- 
ture process zone around the cracks. Additionally, Fig. 
|B^ indicates the existence of anomalous scaling of rough- 
ness as soon as the branching of the cracks is allowed 
[3a |. However, removing the jumps in the crack profiles 
eliminates this anomalous scaling as can be seen from the 
collapse of power spectra shown in Fig. [BJs. Thus it ap- 
pears that the anomalous roughness is two-dimensional 
fracture simulations arises due to crack branching and co- 
alescence of microcracks, which induce jumps in the crack 
profiles. It should also be noted that anomalous scaling 
of crack roughness is readily evident in fracture simula- 
tions obtained using fuse lattice systems. The reason for 
this greater propensity to exhibit anomalous scaling in 
fuse lattice systems appears to be due to scalar nature of 
fuse systems (anti-planar shear model) , which readily al- 
lows for crack branching thereby inducing jumps in crack 
profiles. 
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FIG. 6: (Color online) Scaling of crack widths and power spec- 
tra for the simplified RBM model, (a) Crack profiles (with 
jumps) seem to exhibit anomalous roughness scaling. How- 
ever, the local roughness exponent is same as that noted in 
conventional RBM. The inset shows the scaling of power spec- 
tra, which are parallel to one another but do not collapse onto 
each other, (b) Collapse of power spectra of crack profiles ob- 
tained after removing the jumps in the profiles. Anomalous 
scaling of roughness in this simplified RBM appears to be due 
to these jumps in the profiles. 
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